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n^^TEODUCTIOK 


The  type  of  sandwich  construction  which  is  considered  here  consists 

of  two  thin  anisotropic  Kir chhoff- Love  shells  (facings)  separated  by  a 

three  dimensional  orthotropic  medium  (core)  in  which  the  in-plane 

stresses  are  zero,  see  figure  1.  Since  ^  O  in  the 

core,  only  the  transverse  shear  resultants  S  and  the  mean  normal 
33 

stress  need  be  considered  when  dealing  with  the  core.  On  an 

"  oc 

element  of  a  facing  (see  figure  l)  the  force  and  couple 

-  ^  QC* 

per  unit  of  coordinate  are  evaluated  at  the  surfaces  vzhich 
are  common  to  both  a  facing  and  the  core  (interfaces). 

The  prefix  n  stands  for  O  or  X  according  as  the  quantity  is 
associated  with  the  upper  or  lower  facing,  respectively. 

To  avoid  considering  continuity  of  displacements  at  the  interfaces, 
the  interface  displacements  (^V  )  are  utilized.  In  the  formulation 

of  the  theory,  the  sums  and  differences  of  the  interface  displacements 
(,X>cyp  and  )  are  introduced. 

The  dimensionless  surface  coordinates  are  assumed  to  be  lines 

of  curvature.  Hence,  the  metric  tensors  (  n^oc^  ^ 

coefficients  of  the  second  fundamental  forms  (  and  ) 

are  diagonal  matrices.  Also,  the  coefficients  of  the  second  fundamental 
forms  are  associated  with  the  curvatures  of  the  surface  under  consideration 
The  differential  equations  governing  the  small  deflections  of  the 
above  described  sandwich  shell  are  derived  from  the  Hellinger-Reissner 


Usual  tensor  notation  prevails,  see  reference  1. 
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variational  theorem  2  The  equilibrium  equations  for  the  composite 

shell  are  similar  to  those  obtained  in  3  and  the  boundary  conditions  for 
the  individual  facings  are  comparable  to  those  obtained  for  a  homogeneous 


shell  in 


Ln  4  . 


The  stress  resultant-displacement  relations  for  the 


composite  shell  obtained  here  have  not  been  presented  before.  Repre¬ 
sentative  equations  for  a  sandwich  shell  with  a  viscoelastic  core  are 
displayed. 

The  equations  presented  here  are  applicable  to  plates  as  well  as 
shells 5  however,  they  vrill  not  be  specialized  in  their  general  form 
since  a  complete  theory  of  sandwich  plates  has  been  given  by  G.  A. 
Wempner  and  J,  L.  Baylor  5  • 

Many  authors  have  used  variational  principles  in  the  derivation 


of  sandwich  shell  theories,  E.  Reissner  6  and  C.  T.  Wang  7  ^sed 


the  principle  of  minimum  complementary  energy  to  derive  the  stress 
resultant-displacement  relations  for  a  composite  shell.  Both  Reissner 
and  Wang  regarded  the  facings  as  membranes.  Equations  which  include 
the  bending  stiffness  of  the  individual  facings  have  been  derived  by 
E.  I.  Grigolyuk  8  and  R,  E.  Fulton  9  f^rom  the  principle  of 
stationary  potential  energy.  A  non-variational  derivation  of  sandwich 


shell  theory  is  given  by  Wempner  and  Baylor  j^3J* 

Presented  here  is  a  theory,  developed  from  the  Hellinger-Reissner 
variational  theorem,  which  includes  bending  resistance  and  dissimilarities 
of  the  facings.  The  resulting  equations  are  applied  to  examples 
illustrating  the  effects  of  anisotropic  facings  and  a  viscoelastic  core 
on  sandwich  shell  behavior. 


"'“^Numbers  in  brackets  refer  to  the  bibliography  at  the  end  of  the  report. 


A  THEORY  OF  ANISOTROPIC  VISCOELASTIC  SANDWICH  SHELLS 


1.  Stress  Distribution  thru  the  Core 

In  what  follows  it  is  assumed  that  the  components  of  the  displace¬ 
ment  vector  and  their  derivatives  are  infinitesimals  of  the  first  order 


and  the  squares  and  products  of  these  infinitesimals  are  neglected  when 


compared  with  their  first  pov/ers. 

The  core  is  weak  in  the  sense  that  it  only  resists  transverse  shear 
and  transverse  normal  stresses,  i.e.  =  O.  Upon  setting 


O 


in  the  equilibrium  equations ^  the  core  stresses 


become  statically  determinate 


gives 


Integration  of  the  equilibrium  equations 


and 


[2ul  [j-a:e%zV 


where  CJ  and  are  proportional  to  physical  stress  and 

physical  stress  resultants,  respectively. 

On  the  edge  of  the  core  the  shear  stress  distribution  is  a  priori 
statically  determined  in  terms  of  the  shear  resultant.  From  equilibrium 


4 


of  a  boundary  element,  the  shear  resultant  on  the  edge  of  the  core  is 

This  shear  resultant  must  be  assigned  on  the  edge  of  the  core. 

2.  Core  Stress-Strain  Relations 


The  relative  displacement  of  two  particles  on  the  normal,  one  at 
each  interface,  is 


Presuming  the  core  to  be  orthotropic  with  respect  to  the  surface 
coordinates,  the  stress  and  strain  components  are  related  as  follows; 


5 


Because  of  the  displacement  assumption 


(y) 


and 


'30C 


Substituting,  in  turn,  (?)  and  (5)  into  (3)  and  (7),  (8),  (?)  and 

(6)  into  (4),  expanding  the  integrands  in  ^  power  series  and  neglecting 

2 

^  when  compared  to  one,  there  results 


ad 


31?- 


and 


\  Coc), 


4- 


a' 


6L.E: 


-h 


-h 

3Er  ^ 


JS-LE 


15LE  ^ 


-i- 


6 


-i(8hh,«-ie,«)].  C^°) 


Equations  (9)  and  (lO)  are  the  core  stress-strain  relations. 

3.  stress  Distribution  thru  a  Facing 

The  force  and  couple,  per  unit  of  coordinate,  on  an  element  of  a 


facing  are  (see  figure  l) 


OC 


“■|®7=r9 


vhere  Lj,m°^ood  are 

physical  stress  resultants. 


proportiona.1  to 


Neglecting  terms 


of  order  ^  0^  y  stress  resu.ltants  are  related 


to  the  stresses  as  follows; 


a  c,  jX°'^  ^  ^  ^ 


Z —  n 


93 


((13) 
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notice  that  and  ere  sy.mrnclric . 

Guided  by  (ll)  and  (.12),  the  stress  components  are  presumed  to  have 
the  following  form 


3_ 


-h  2^0^411 


6^rn'^ 

/  4  oa 


(lA) 


or 


(is) 


The  normal  stress  is  assumed  to  be  zero. 


4.  Strain-Displacement  Relations  for  a  Facing 

Tlie  facings  are  presumed  to  be  thin  Kir chhoff- Love  shells,  i.e. 
normals  remain  straight  and  normal  to  the  interface  surfaces. 

If  extension  of  the  normal  is  neglected,  the  displacement  of  a 
particle  in  a  facing  is 


^  ~  ry  (n^3~ 


The  deformed  and  undeformed  unit  normal  vectors  are  related  as 


follows 


5 


8 


A  —  A  -I-  — 

-  Q3  +  nO.  . 


(17) 


Substituting  (l?)  into  (16),  one  finds 


V  ==  -hnO^ 


The  covariant  components  of  the  displacement  vector  are 

w  =  V 


(J6) 


^o(.~^rP  a^o<3 


V,  -  3.L.  ± 


n  vj 


(iq) 


(20) 


=  L  (1  +  ±'^)  (2i) 


^3^  =  ).  (22) 


Equations  (20)  and  (2l)  are  obtained  directly  from  the  definitions  of 


X>t/p.  and 


Because  of  the  displacement  assumption  1 


SL'^ccU-'-  \slo< 


(33) 


uses  (18),  (19),  (20),  (21)  and  (22)  in  (23)  and  neglects  terras 


of  order  n^l  j  then 


9 


= 


X"  _  '~)X 

r>  n  <=><{3 


(24) 


where 


Y  -  X 

CL  2 


p^j/3  G 

-  2  ^^aVa  1  Css) 


and 


>{  =-^- 
D  £ 


rfS!3<^^'^o‘^3iejx~  (26) 


Equations  (2^)^  (25)  and  (26)  are  the  same  strain-displacement  relations 


obtained  in 

o 


n^o<& 


with  Q  and  yO 

o<3  y^> 


of 


10 


and  YK  ..  as  defined  hers  agree 

U.  o<f5 


<=^  y  ocf3 

The  shear  strain  is  zero  at  the  interfaces  and  will  he 

3oC 

assumed  zero  throughout  a  facing. 


5o  Hellinger-Reissner  Three  Dimensional  Variational  Theorem 

The  equilibrium  equations,  stress  re suit ant -displacement  relations 
and  boundary  conditions  for  a  sandwich  shell  will  be  derived  from  the 


'ollowing  variational  principle  of  Hellinger  and  Reissner 


2  . 


The  state  of  stress  and  displacement  which  satisfies  the  differ¬ 
ential  equations  of  equilibriinn  and  the  stress  displacement  relations 
in  the  interior  of  the  body,  and  the  conditions  of  prescribed  stress  on 
the  part  and  of  prescribed  displacements  on  the  part of  the 
surface  of  the  body,  is  determined  by  the  variational  equation 


10 


6 .  Contribution  to  the  Hellinfier-Reissner  Theorem  from  the  Facings 


The  normal  stress  and  the  shear  strains  ^ 


have  been 


assumed  zero.  Having  zero  shear  strains ^  j  v/hile  there  exist  non¬ 
zero  inpiane  stresses  ^  and  non- zero  shear  stresses 

requires  the  elastic  coef f icients^^_2cx^ 

zerCo  Thus  the  only  contributions  from  the  facings  to  the  volume 


integral  of  the  variational  theorem  are 


Substituting  (15)  and  (24)  into  the  integrals  thru  a  facing  thickness 
and  neglecting  terms  of  order  ^  CX  ^  one  finds 


j  ^3_  dn°'^±BLn^  y 

-J 


11 


and 


(30) 


drn'^^  ±2Lm°'^ 


2L2d„d2  J 


d'-'  n 


n 


1 


j^n°^±  SLn^JQn  ^9±  2  L  r/?J  -f 


+ 


4- 


2lF3%3^  pn  °^±  2Ln'^ 


''P=t2Lm^9j  -f 


4- 


2L^J%d^^j 
3 


|dr)^ir2Lr-)^7 

nL  J 


-i 


J  rrf^^+i2L rrf^  -h 


dm‘^±2Lm’^ 


pm  ^7±  S-Ltr/ll^si) 


The  sinrface  of  a  facing  consists  of  three  parts;  the  interface 
surface,  the  exterior  face  and  the  edge.  Over  the  exterior  face  and 


the  interface 


—  Oj  i.e., 


stresses  are  prescribed.  The  integrals 


p  vp.  J^S , 


n 


over  the  interfaces,  from  the  facings  are  the  negatives  of  the  corre¬ 
sponding  integrals  from  the  core,  consequently  they  sum  to  zero. 

The  load  on  the  exterior  face  of  a  facing  is 


-  [nP 


r’d^oc+a  P  °-s  J]j 


czJ 


(32) 
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where  and  are  proportional  to  physical  force 

per  unit  undeformed  area. 

Using  (18),  (19),  (20),  (21)  and  (32),  neglecting  surface  load  times 
rotation  terms  and  neglecting  terms  of  order  and  one 

obtains 


J'l'pV^c^s  = 

S’  S  ^ 

D  ^ 


oc  oc 


From  (13),  (l^i),  (l8)  and  (19),  after  neglecting  terms  of  order 
and  it  is  found  that 

}[pV,4n+JT(v,-v,)p'‘4^  = 

_o . 


D  2. 


^  rf—  9=i“^\/  —  —  rn°^V  ■+ —  'q^V  — 


n^i 


d-  r  w  ,  * 


c  L 
£?  ^2 

“  rC  f-' 


2  oc  /  .  .  \  nLioc  J  y) 

d-^nQ  (oVs-a'^s  Jj-^  ■ 


(34) 


This  is  the  contribution  to  the  surface  integrals  from  the  edge  of  a 
facingo 

7o  Contribution  to  the  Hellinger-Reissner  Theorem  from  the  Core 
For  the  core  the  equilibrium  equations  have  been  identically 
satisfied,  the  stress-strain  relations  have  been  determined,  the  boundary 
condition  has  been  obtained  and  over  the  interface  surfaces  has 

been  presumed  zerOo  Thus  the  only  contribution  to  the  variational  theorem 
from  the  core  is 


The  outward  unit  normals  to  the  core  interface  surfaces  are 

Qn  =  ±aL^'^  (35) 

Equations  (l),  (2),  (20),  (21)  and  (35)  give 

.(6^±£^y^  (±i+afct^)(5^± 

Js.  (36) 

8.  Hellinger-Reissner  Variational  Theorem  for  a  Sandwich  Shell 

Upon  substituting  (28),  (29),  (30),  (31),  (33),,  (34)  and  (36)  into 
the  variational  theorem  (27)  and  using  Green's  theorem  1  ,  one  obtains 
the  following  variational  equation  appropriate  for  a  sandwich  shell. 
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5 


fVr)  \  i_ 

-,  Wn^9-  -  j- 


15 


-'■’-  f ^'i’- 


<>9)1+5^15 


-p‘+xy] 


-h  Syixrx 
oC 


-T^  n^) + 


16 


c^3  aVs 


17 


In  the  boundary  integrals  of  (37)  the  following  moment 


equilibrium  equation  was  used; 

This  equation  is  derived  in  |^3  • 

Equation  (37)  is  the  required  variational  equation  for  a  composite 
sandv;ich  shell. 


9o  Equilibri^Jin  Equations 

The  Euler  equations  resulting  from  operating  on  (37)  and  corre” 
spending  to  and  are 

+  ^('T^rn^+Xi^"^)llr  =  0  C38) 

and 

+  =  o-  (39) 

Forming  linear  combinations  of  the  Euler  equations  associated 
with  and  one  finds 

^  "  r/;«  -  a  "  +  a  T '  B  ‘X  + 

-am«'/i:V;b"+Gf>;('t")']-o  fao) 
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and 


- S  -  +  n n^<^- 

-a-»rnWfrJ'fa"+-r;('fa2fJ-o.  (41) 

Equations  (38)?  (39)?  (^*3)  and  (4l)  are  the  equilihrium  equations  for 
a  composite  sandwich  shell.  If  the  equilibrium  equations  of  3  sre 
specialized  to  small  rotations  and  ii  is  neglected  when  compared 

to  one,  the  resulting  equations  are  the  same  as  (38),  (39)?  (^=0) 

(41). 

Equations  (38)  and  (40)  are  identified  with  the  equilibrium  of  a 
gross  element  of  the  composite  shell. 

10.  Stress  Re suit ant -Displacement  Relations 

Combining  the  Euler  equations  corresponding  to  S  H  S  ^ 

Srn^^  and  in  a  suitable  way,  it  can  be  verified  that 

+  Ol(oc1^-O^T^  --^4^!^)]  + 


19 


Co<)  3  J 


Ji  [-CC -^tU^)jl^  + 

-f  5c  (a^  -:xh%'fiZ^)jy+  “^C/  ^ 

■'■  )('^  -^4Ij  ^) 

'•■  i  2  Ji  J  2^  O^jOi  ■'■ 

+-^^t\)ll^  -  ^ 0^3JX+-^I=>%]\  + 


oc^  '  '--r  -^oc^  JK^''3j}i'  '  y^C^)  , 

^[l\  )J 

Jl  J  O ^ 


-7^ 


OC 


6x')\ 


‘Co^J  y//^ 


3j^ 


/  (cx)  \ 

■^^c^))y-^ 


y^cyj 

i(d) 

y^{eri. 


,  C-^s) 


n^'^=2  JiJo^X  A-:X  B  f - oc(m^  - 


20 


21 


d)//^  +■ 

-^a^T^-o^  bJ5  jj,  (43) 

abg;^)] + 

+  kajLjA'^l-^  (z^  - 

-^^2'^)//^-/S(^-^^t;^)|^  - 
“/^4^ 

+ JiJ  -h 


22 


m 


9 


23 


3  X  L'’  B-^'^°‘^jj!' -+- 

+a(-y^+^F'r;fX^.^+^^-b^S)  + 

s  a  Ji.  i- 


24 


Remembering  that  and  have  been  neglected  vhen 


compared 


to  one,  care  must  be  taken  vhen  using  ('12),  (-G),  (44)  and  (45 )  since 

and  ^  are  contained  in  OC  ^OC ^  ^  ^  ^  ^ 

may  be  contained  in  J3  ^  . 

Terms  multiplied  by  OC^  ^  and  ^  when  =  ^  by 

V7hen  the  facings  have  the  same  physical  properties,  by 

OlIdT\  and  by  Qf  ^  and  DfV^  are  due  to  the  variation  in 
(p^j  ^j3 

the  geometry  thru  the  composite  shell  thickness*  If  the  sandwich  snelL 
is  thin  these  terms  can  be  neglected 0  Hence  for  a  thin  sandwich  shell 
V7ith  equal  facings,  (42)  reduces  to 

-  i  L?  If  4<  I  )//^  ] . 

1 

This  is  the  same  stress  re  suit  ant -displacement  relation  obtained  in 

Equations  (9)3  (lO)?  (38),  (39)>  (4o),  (41),  (42),  (43),  (44)  and 
(45)  form  a  system  of  21  simultaneous  differential  equations  in  the  -rr- 

v&rx&hies  s°^  rh'^i^oncl 

11.  Boundary  Conditions 

The  boundary  condition  for  the  edge  of  the  core  has  already  been 
given,  i.e.  the  shear  resultant  on  the  edge  of  the  core, 


^ 

must  be  specified. 

Since  on  a  normal  to  the  core  mid- surface  at  the  edge  of  the 
composite  shell  stresses  may  be  prescribed  for  one  facing  v.hile  dis¬ 
placements  are  prescribed  for  the  other  facing,  the  boundary  conditions 
for  the  individual  facings  will  be  given.  Using  1  3] 
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integrating  by  parts  and  then  setting  the  resulting  coefficients  of  the 
varied  quantities  in  the  line  integrals  equal  to  zero,  the  required 
facing  boundary  conditions  are 


m  iJ  LU  =  py^  I J  I J 


n 


Oi-6) 


and 


^o</S~l  I  nU-oc 

■  nj  -  J  nj 

r u  ^c<  cx 

L  ^  .  I  _j  £2  (J  ”  ^ 


£?  ^  O  U 

on  ^C.  and 

D.  ~  o  J 


n  Vs 

and 

~  n  ^3j>C 

These  boundary  conditions  have  the  same  form  as  those  obtained  in 
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12.  Stress-Strain  Relations  for  a  Sandv.ach  Shell  V7ith  a  ViscoelaGtic 


In  the  following  two  sections  a  sandwich  shell  with  a  viscoelastic 
core  is  considered «  Representative  equations  for  this  shell  sre 
displayed. 

Only  the  core  v/ill  be  presumed  viscoelastic,  however,  viscoelastic 
facings  could  be  treated  in  the  same  way. 

The  core  stress- strain  relations  are  altered  as  follows ;  ^  and 
'  are  replaced  by  /y  El  and  EL  ^  respectively,  and  all 
:ther  functions  of  time  are  replaced  by  their  Laplace  transforms,  e.g 


(9)  becomes 

^  1^--^  L — 

^ 


Sine, 


[v"+ 


are  nvot  functions  of  time,  the  stress 


resultant-displacement  relations  for  the  composite  shell  are  converted 
simply  by  substituting  Laplace  transforms  for  all  time  functions.  To 
illustrate  this  a  few  terms  of  the  equation  corresponding  to  (42)  are 
presented; 
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-f - . 


13 »  Equilibrium  Equations  and  Boundary  Conditions  for  a  Sandwich  Shell 
with  a  Viscoelastic  Core 

The  equilibrium  equations  for  the  composite  shell  and  the  boundary 
conditions  for  the  individual  facings  and  the  core  are  obtained  by 
merely  replacing  all  functions  of  time  by  their  Laplace  transforms 3  e.g. 
(38)  and  (46)  beccme 

and 

Li  LX  =■  ^/^^LX  LX 

It  has  been  assumed  that  arid ^  ^ oc  independent  of 

time , 

EXAMPLES 

The  theory  presented  here  is  valid  for  sandwich  shells  (plates) 
with  thin  Kir chhoff- Love  shell  (plate)  facings.  However,  the  following 
three  examples  are  only  concerned  with  sandwich  shells  (plates)  with 
membrane  facings.  The  facings  are  presumed  membranes  so  that  the 
influence  of  a  hole,  of  orthotropic  facings  and  of  a  viscoelastic  core 
on  the  behavior  of  a  sandwich  shell  (plate)  can  be  studied  without 
unduly  complicating  the  examples. 
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.l4.  Circular  Plate  v;ith  a  Circular  Hole  at  the  Center 


Consider  a  simply  supported  circular  plate  with  a  circular  hole  at 
the  center  loaded  by  a  uniformly  distributed  bending  couple  around  the 
outer  boundary.  The  facings  are  isotropic  membranes  with  similar 
physical  properties  and  equal  thichnesses.  The  core  is  presumed 
isotropic.  The  dimensionless  surface  coordinates  are 

where  p  and  (j)  are  polar  coordinates  (see  figure  2). 

Dae  to  the  symmetry  of  the  plate  and  the  applied  edge  couple, 

and  the  remaining  dependent  variab.les  are  independent  of  0  . 

The  equilibrium  equations  which  are  not  identically  satisfied  are 


0 

II 

to 

n 

.  ~h 

=  O 


11 


(4-1) 

C4-S) 

Cso) 
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The  cere 


?  tr e s  u rain  re .l.a t i on s  ar e 

aL_  v33 

EE  ^  ^ 


,+  2i_q 


CL^ 

&LE^  I 


(si) 

Cse) 


Th^^  streiss  re  suit  ant- displacement;  relations  for  the  composite  plate  are 


^22.  _  s 


CLR  ^e:  u  /  . 


'c-{'AA-'0~^  1 -/— 


(-53) 


r)^‘=a^a/RT^£ 


2.5*- 


6^-^^  J  1  ■'■  ■!  )J. 


The  Doundaiy  conditions  are 

[" ‘i‘-, = 


3],  (■**) 

(57) 


(5S) 


(60) 


0^1  "=j 


=  0. 
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From  (48)  and  (5I)  one  sees  that 


o. 


The  solution  of  (47)  and  (57)  is 

<5  ^  =  O. 

Substituting  (53)s  (5^)  and  (6l)  into  (49)  gives 

1  /-i 


=  o. 


Equation  (63)  and  the  boundary  conditions  (58)  yield 


i 


=  o 


(6i) 

(^62) 


(63) 


(64) 


Hence 3 


O 


In  the  same  way  (64)  was  obtained,  from  (50),  (62),  (52),  (55)^ 
(56)5  (59)  and  (60)  one  finds 


yvi 


a  -  r r 

This  solution  has  exactly  the  same  character  as  the  solution  of  a 


homogeneous  plate 


12 


If  the  facings  had  been  thin  plates  instead  of 
membranes  the  problem  would  have  been  greatly  complicated  and  the 
character  of  the  solution  would  have  been  different.  The  character  of 

the  solution  would  depend  on  the  boundary  conditions,  however, 

—  I  _  ^ 

and  S  in  general  would  not  be  zero  and  AM-  v/ould  be  considerably 

more  complicated, 

A  sandwich  plate  with  equal  facings  and  a  hole  (circular  or  not) 


32 


loaded  by  the  same  inplane  edge  tensions  on  each  facing  has  exactly  the 
same  solution  as  a  homogeneous  plate.  In  this  case  the  facings  can  be 
either  membranes  or  thin  plates. 


=  OC.  s 

= 

I 

£ 

i^'j 

2^1 

II 

> 

1^2. 

o, 

~  /^a  i 
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FIG.  3,  SQUARE  PLATE 
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(65) 


Since  the  plate  is  simply  supported  and  symmetric 

O  1 


n  =  o  j. 

The  equilibrium  equations  not  identically  satisfied  are 

L2^33  p  =  o_, 

O- 

The  core  stress- strain  relations  are 

“  iOO "  ^2Z_G  6L-B 

The  stress  re suit ant- displacement  relations  for  the  composite 
plate  are 

-4W  +  S  )]  , 


(66) 

(67) 

(68) 


(69) 

(70) 


(7J) 


(72) 


The  boundary  conditions  are 


(Oi.oc'j 


.o,-°'We'-o.rO-  72) 
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From  (67)  it  is  seen  that 


is  a  constant,  hence,  from  (69) 


it  follows  that  is  a  constant.  According  to  (66),  'S 

is  a  constant  so  that  (70)  reduces  to 


I  __S^ - 

ioo  ^  SL-G 

Using  (65)5  (71)  and  (72)  in  (68)  one  obtains 

4-1  fpriJli-— .  =J122-_  ~1  _ 

4-ooo\_ 

_  ^  J-  )  ■^■^^22.  "1  I 

SLGU 

-t-  ■*  =  O 

4-i  r=22.n  =222a “* 

^ooor  ^ 


674.) 


(7^) 


SLGr  Jje  ^  ^ 


Jaa 


+  =  O. 


676) 


where 


^(ccocX^^)  _  ^  C<Xo<X^f3) 


-B 


-B 


[g  Jaiig222a_(^g  jjssj)® 

<ro».)2sr6 ^ixaSgUO^)  ■ 

gijii^22ea__^^jjaaj^a 
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Equations  (66),  (75)  and  (76)  are  three  simultaneous  differential 


equations  in  the  three  dependent  variables 


S  ctnd  -S 


The  following  series  satisfy  the  boundary  conditions 


L.CI 


■■  ^  Ak  sin(r^0^^sir)(snr9^)  . 

(s'rro'^)  ^ 

r4s=ocld 


S  ^ 


where 


Q  = 


S  C  sin 

lr^S=:  odd 


P 


(r^yG  ^)cos  (s'rr 


Substituting  these  series  into  (66),  (75)  and  (76)  one  obtains  three 
simultaneous  algebraic  equations  in  the  three  sets  of  constants  ^  ^ 

and  C_  ,  Solving  these  equations  one  finds 

ris  r~s 
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II 

VO 

< 
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The  remaining  unknow  functions  can  now  be  determined „  One  obtains 
from  (7^)  and  from  (68)  and  (73); 


l:q 


^  rTfA^^cos(r'iYe^)stn(s7T'&^) 

^^r4s  =  oJcJ 


c.os(r7YG^'^ s/n  (s'TYG'^ 


r-is=zaei<i 


L-Cl 


lOO 


Z.  sTr'A^^sin  ^r^-Q^'^cos(s7r'0^^  -h 


odd 


so  S*  sin(r~^0-^\cos~(s'Tf^'^ 

ris  =  odd  <-s  ^  J  K  y 


(t^G^^s]n  (s'frO’^ 


Figures  4^  5j  6  and  7  show  disp-lacements  and  stress  resultants  for 
a  sandv/ich  plate  whose  upper  facing  remains  isotropic  while  its  lower 
facing  ranges  over  various  degrees  of  orthotropyo  equal  to 

Young’s  modulus  for  the  upper  facing  and  .E^  decreases  from  . 

The  cross  modulus  is  assumed  to  have  the  following  form 


i-^iS  V  J  ^2.  • 
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From  the  figures  we  see  that  as  the  lower  facing  ranges  over  various 
degrees  of  orthotropy  all  displacements  and  stress  resultants  hehave  as 
one  would  expect.  The  stress  resultants  S  snd  are  larger  than 

and  ID  since  the  plate  stiffness  in  the  direction  is 

o 

greater  than  the  stiffness  in  the  ©  direction.  For  the  same  reason 
the  rotation  is  less  than  • 
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i-Q,  Vg. 


rr*o 


0 


4 


FIG.7,  BENDING  MOMENTS 


l6.  Infinite  Circular  Cylinder  with  a  Viscoelastic  Core 

In  order  to  study  the  effect  of  a  viscoelastic  core^  an  infinite 
circular  cylinder  loaded  by  a  concentrated  uniform  ring  load  acting  at 
©  is  investigated.  The  facings  are  isotropic  membranes  vith 

the  same  thickness  and  physical  properties.  The  core  is  isotropic 
with  an  infinite  Young *s  modulus  in  transverse  extension. 

The  dimensionless  surface  coordinates  are 

=  ©^= 

see  figure  8.  Hie  assumed  viscoelastic  character  of  the  core  is  that 


of  a  standard  linear  solid  as  shown  in  figure  9,  For  this  example 


we  take 


1 

so  ^ 


1 

1600 


In  the  sequel  the  variation  of  the  geometry  thru  the  thickness 
of  the  composite  shell  has  been  neglected. 

From  s^-mmetry  of  the  shell  and  the  load 


and  all  remaining  functions  are  independent  of 
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FIG. a,  CIRCULAR  CYLINDER 


49 


The  following  equations  are  the  time  Laplace  transforms  of  the 
equilibrium  equations,  the  core  stress- strain  relations  and  the  stress 
resultant-displacement  relations  for  the  composite  shell 


-^P6{9V  O  , 

(77) 

-*n^'=0, 

IN 

N^ 

*n%-0^ 

(79) 

*s^  —  0_, 

(so) 

—  SEJ  '  ZLrfQ 

(SJ) 

728  , 

(82) 

^  J  J  -  Qi^  Fo’^y  /  _ /T^ 

^  iJ6480oL^^^j^  1600  3^ 

1J648oL  J600  3^ 

aa],  (S3) 

Notice  that 


JO 
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re  J  1 


s  the  Fourier  trar;f?fGrin  parameter  and  a  \/  over  a  function 


indicates  a  Fourier  transform* 

The  inverse  Laplace  transform  of  (89)  is 

r  5?^-/- 
\  1 64- 


_  _  128000 

LQ  “ 


^4-  /  1000  I  a3a9feO 

j  -6-  4  ^ 


r  5-  e2  ,  J4^56 

J64-^  -o  -5:] 

-  — » 

e4_f_  200  e^_i_  23296-0 


^4.  iQQO^g,  e3e960 


S'Xp 


232960 


'/V  ^^VT) 


("90) 


L2.-E 


^  = 
J  ^ 


h serving  thau  yUJ-X  is  an  even 


1  function  ^ 


it  IS  seen  liiat 


AAJ~^  _  S. 

us-y^ 


oo  y/ 

La 


{g)cos(j&^)Jj  (m) 


;xGendinB  the  exponential  function 


Q-ycp 


JOOOroa,  238860 

3  ~7rr~  s  ^  j^i 


-t 


4-00'^!  ^4,  200^2,  232960 

VJ  :4l 


(92) 
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in  a  power  series  and  comparing  the  integrals  which  result  from  substi¬ 
tuting  (90)  and  the  power  series  of  (92)  into  (91)^  it  is  seen  that  in 
approximating  (92)  by  exp  C  ^00^  only  a  term  of  order  icr^  as 
compared  to  one  is  being  neglected  vrhen  4  12.00,  With  this 

approximation  (91)  reduces  to 

^  -  3oao(^i_ 

•  [cos  (5. 6aa^)  -/-  O.  feos  s/'n  (5.  629^)]  - 

-^0.907s;n(f6..04-&^)J.  (93) 

Making  the  same  approximation  in  the  integral  form  of  and 

satisfying  (85)  one  finds 

=  J20— i2o4— 


lcos(s.62.0^^  —  0.164  sin  (s.  62©®)J  — 

_  i20e“'i^e  - 

-  0.030.4  sin  (6.040^)1. 


(94) 


The  stress  resultant  /O  can  he  determined  from  (87)  and  (93) 

^  a 

To  determine  3  one  uses  (77),  (8?)  and  (86)  to  obtain 


^  =  S(e^)  -  J- 
P  ^  240\.  L.Q.  J' 

Integrating  (95)  snd  using  the  boundary  conditions  on  and 

one  finds 


C9S) 
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Equatica  !3l)  can  b?  written  _ ^ 

+  (’7) 

Again  aiproximtlag  i92}  by  exp  beccmea 

£  ^  _  ir 

^  =  [4-.  fe&  (1  -  e^)  -loo 
•  e  ^•^®®Ls(5.62.&^)  -  Ifse  (1-  e^)  - 


eoi.  i_a 

I 


+|Q623-t-iOo('^5^)  e  ^•'^5s  (4.04-0^)- 

-  [377+3.24.(1^;^)  e^^e^'^sln{6.0^&^).  (f!^ 

Equations  (83),  (8l)s  (96)  and  (86)  yield 

#= i]i+8o  5  ^  ")  - 

4i  A^,aa\  _jL  jL./(:d±^\  J— 


'1600V  L.GI 
,  /-t 


siLa 


_l_  — - 
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it  i 


In  evaluating  ^ 

3j  c.'S. 


s  seen  that  it  contains  and 

that  (92)  can  no  longer  be  approximated  by  exp  ^  ^^00''^  Here  one 
must  use 


eocp 


-t 


1000^^2. ,  '232.^ (=>0 

j  -znr^j  — 


4007:: 


—  ^ 


J- 


'7f 


Boo  .  2.32960 

pp 

-I 


After  some  manipulation  one  finds 


JOOl 


'n 


^  -  [s.ya  (i  -  © 

•jf  -n  y  _ _ /P. 


2  ^(400 '^r)  e  os  (^^r.  6  2  -f- 


4- 


^Asfi-e  4-  o."7 vs (43^) s  ■^-oo'^n. 


4 


j 


•  e  ^'^'^%!n(s.i>'S.&^)  -  [4-28  +  J.10^~^)] 

.  +  [4.42  - 


.00/; 


4-007^)  ^ 


-6 _ 

4oo^ 


~  6-2^©^ 


in  (99 


From  (78)  and  (84)  we  have 

-  5^0^)  +  24i-S-^V  ^iOO\ 


Equations  (93),  (9^)>  (98)  and  (99)  are  only  valid  for  ^  12.00 

and  For  it  is  observed  that  and  n  ^2. 

are  even  functions  of  ^  ^  and  that  JjU-^  and  si’s  odd  functions 

ot  Equations  (96),  (87)  and  (lOO)  are  valid  for  all  zf  and 

For  zf  ■=  00 


=  —  30£0e  jcos(s.b2e^^  -h 

“f“  0. 60£  s/o  6>£©  ^3  ^  j 

^  =  i20  —  S20e^'^^^^\cos{s.(:.Ze‘'^)  — 

—  O.  s/n  {^.££©^3^  j  (iOj3 

=  4-.66  e"'^'^®®Y^os(;^.6£©®3  - 

—  ss. o sir-, (s.b.s.e^)J  j,  (jos) 

-  3.V3  e  ^■^®®^|cos(^5'.6£©®)  - 

-  J.  13  sin  {S.6Z&^)1. 


=  ISO  —  iSO  e 


^  =  4-66  e  /ccps 


<-rT  r— - 6.^2©* 

100(J^  )  —  -  e 


Equations  (98)  and  (102)  show  that  is  discontinuous  at  ^  ^=:  O 

which  does  not  agree  with  physical  reality „  From  (83)  it  is  seen  that 


>OU±  must  be  discontinuous  if  n 


to  finite  at  ^ 


If  the  facings  had  been  thin  shells  instead  of  membranes  this  incon¬ 
sistency  “would  not  have  arisen^ 

As  can  be  seen  from  (9^)  snd  (lOl),  at  is 

independent  of  . 
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In  figures  10  to  only  the  factions  for  t  equal  zero  and 
infinity  are  plotted =  The  functions  at  ^  =  1200  are  so  close  to  the 
functions  at  “(t  =sOO  that  they  are  almost  indistinguishable  on  the 
figures. 

For  the  nmerical  values  of  the  physical  constants  chosen  the 
displacements  and  stress  resultants  do  not  vary  greatly  as  functions  of 
time.  However^  for  a  different  set  of  numerical  values  for  the  physical 
constants  this  may  not  he  true. 
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FIG.  11.  GROSS  AXIAL  DISPLACEMET 


NOTATIONS 


The  tensor  notations  of  j^2j  are  utilized,  Latin  suffixes  take  on 
numbers  2  and  3  while  Greek  suffixes  take  on  numbers  1  and  2, 

Repeated  indices  are  not  summed  when  enclosed  by  parentheses.  The 
prefix  n  stands  for  O  or  X  according  as  the  quantity  is  associated 


with  the  upper  or  lower  facing,  respectively.  If  two  signs  appear,  i^e. 

-f- 

—  ^  j  upper  (or  lower)  sign  applies  whenever  reference  is  being 
made  to  the  upper  (or  lower)  facing,  A  comma  denotes  partial  differen¬ 
tiation,  i,e,  ^  vertical  bar  (|)  denotes  covariant 

diiferentiation  with  respect  to  the  three  dimensional  space  while  a 
double  vertical  bar  (  j| )  denotes  covariant  differentiation  with  respect 


to  the  core  mid- surface  coordinates. 


Symbol 


Description 


a  characteristic  length  of  the  core 


mid- surface 

thickness  of  the  core 

thickness  of  a  facing,  O  —  O  or  I 

dimensionless  siirface  coordinates, 


lines  of  curvature 
dimensionless  normal  coordinates 
base  vectors  P.cx  where  p  is 
the  dimensionless  position  vector  of 


the  core  mid- surface 

unit  normal  to  core  mid- surface 
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D.  ^^^3 
1^0 


o.9s' 

a.  — 


£ 

CL^ot^ 

0.^3 


dimensionless  interface  base  vectors 


(see  O,^  ) 

coefficients  of  the  second  fundamental 


form  for  the  core  mid- surface 


coefficients  of  the  second  fundamental 


form  for  the  interface  surfaces 


base  vectors 


^/S 


determinant 


^CK'  ^/9 

determinant  I^CK/S  I 

determinant  /  O  ^1 

permutation  symbol 

deformed  unit  normal  to  interface 

sur^faces 

a_.a  unit  normal  to  core  edge 

o< 

at  the  mid- surface 


nM»< 


unit  normal  to  the  edge 


of  a  facing  at  the  interface 

unit  tangent  to  the  edge 

u.  2 

of  a  facing  at  the  interface 

mean  curvature  of  core  mid- surface 

Gaussian  curvature  of  core  raid- surface 

^Ljsd 
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r-j 

n 


o  '^yn  ^  _ 

sa'T'f 

cx^ 


yihjad^  y?jJj.cl^ 
yLij^d^- 
VL^j^d^  + 
VL^o^d'^-  yuj^d^ 

Christoff el  symbols  of  the  second 
kind  evaluated  at  the  interface 
surfaces 

[C/]e®=o 

ry^  y  ,  ry\  S' 

3  ~*~  i  -^Oc^ 

ry^y  ^  ry^  y 

o  ’^ocfSi 

volume  of  a  body 

surface  of  a  body 

core  mi.d~sixrface 

interface  surfaces 

exterior  faces  of  the  facings 

part  of  CJ^  on  which  the  stresses 


are  prescribed 


£3  J 

a^z 

aC,- 


part  of  c/  on  which  the  displace¬ 
ments  are  prescribed 

for  the  edges  of  the  facings 
<=4  for  the  edges  of  the  facings 


part  of  interface  boundary  curves  on 


which  the  stress  resultants  are 


prescribed 
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1C  c  IE  5  iQj 


C 

Q^'2. 


2. 

2UJ. 


^rs  - 
r^rs 


^33 


part  of  interface  boundary  curves  on 
which  the  displacements  are  prescribed 
dimensionless  arc  length  along  inter¬ 
face  boundary  curves 
dimensionless  arc  length  along  the 
normals  to  the  edges  of  the  facings 


at  the  interfaces 


^  '  j  displacement  vector 

'^r)s  Vs/r* 

interface  displacement  vector 
a  r'  (^tiy )  average  displacement 


of  the  interfaces 

of  the  interfaces 


strain  tensor 


stress  tensor 


relative  displacement 


CLL 

tOc 


L?  r  _ _  33  , 

OL(on'^^  — 

OL(^nr)^^— 
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extension 


the  core 


01 

shear  modulus  of  an  ortho tropic  core 
shear  modulus  of  an  isotropic  core 
Young's  modulus  of  an  isotropic  facing 
Poisson's  ratio  for  both  facings 
elastic  coefficients  for  a  facing 


elastic  coefficients  defined  by 


n 


O  ^  J.  Yr^ 


■kc 


rsk.JL 


hjt 


a  symbol  placed  over  a  quantity 
indicating  that  the  quantity  is 
prescribed  on  the  edge  of  a  facing 
a  symbol  placed  over  a  function  of 
time  to  indicate  the  Laplace  trans¬ 
form  of  the  function 
Laplace  transform  parameter 

=  l+S.-Xh+3.^k. 
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